Fast quantum gates for neutral atoms 
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We propose several schemes for implementing a fast two- 
qubit quantum gate for neutral atoms with the gate operation 
time much faster than the time scales associated with the ex- 
ternal motion of the atoms in the trapping potential. In our 
example, the large interaction energy required to perform fast 
gate operations is provided by the dipole-dipole interaction of 
atoms excited to low-lying Rydberg states in constant elec- 
tric fields. A detailed analysis of imperfections of the gate 
operation is given. 

PACS: 03.67.-a, 32.80.Pj, 03.67.Lx, 32.80.Rm 



In recent years, numerous proposals to build quantum 
information processors have been made . Due to their 
exceptional ability of quantum control and long coher- 
ence times, quantum optical systems such as trapped ions 
1^ and atoms |^ , and cavity QED , have taken a lead- 
ing role in implementing quantum logic in the laboratory. 
Quantum computing with neutral atoms |5| seems partic- 
ularly attractive in view of very long coherence times 
of the internal atomic states and well-developed tech- 
niques for cooling and trapping atoms in optical lattices, 
far off- resonance light traps and magnetic microtraps . 
Preparation and rotations of single qubits associated with 
long-lived internal states can be performed by addressing 
individual atoms with laser pulses. A central issue is to 
design fast two-qubit gates. 

First of all, it is difficult to identify a strong and con- 
trollable two-body interaction for neutral atoms, which 
is required to design a gate. Furthermore, the strength of 
two-body interactions does not necessarily translate into 
a useful fast quantum gate: large interactions are usually 
associated with strong mechanical forces on the trapped 
atoms. Thus, internal states of the trapped atoms (the 
qubits) may become entangled with the motional degrees 
of freedom during the gate, resulting effectively in an ad- 
ditional source of decoherence. This leads to the typical 
requirement that the process is adiabatic on the time scale 
of the oscillation period of the trapped atoms in order to 
avoid entanglement with motional states. As a result, 
extremely tight traps and low temperatures are required. 

In the present Letter we propose a fast phase gate for 



neutral trapped atoms, corresponding to a truth table 
Id) ® |e2) e"i<^2¥'|g^^ ^ ig^) for the logical states je^) 
with 6; = 0, 1, which (i) exploits the very large interac- 
tions of permanent dipole moments of laser excited Ryd- 
berg states in a constant electric field to entangle atoms, 
while (ii) allowing gate operation times set by the time 
scale of the laser excitation or the two particle interac- 
tion energy, which can be significantly shorter than the 
trap period. Among the attractive features of the gate 
are the insensitivity to the temperature of the atoms and 
to the variations in atom-atom separation. 

Rydberg states ^ of a hydrogen atom within a given 
manifold of a fixed principal quantum number n are de- 
generate. This degeneracy is removed by applying a 
constant electric field £ along the z-axis (linear Stark 
effect). For electric fields below the Ingris- Teller limit 
the mixing of adjacent n-manifolds can be neglected, 
and the energy levels are split according to AEnqm — 
SnqeaoS /2 with parabolic and magnetic quantum num- 
bers q — n— 1 — |m|,n — 3— |m|,...,— (?! — 1 — |to|) 
and m, respectively, e the electron charge, and ao the 
Bohr radius. These Stark states have permanent dipole 
moments fi = ^l^Gz = 3nqeaQez/2. In alkali atoms the 
s and p-states are shifted relative to the higher angu- 
lar momentum states due to their quantum defects, and 
the Stark maps of the m = and m = 1 manifolds are 
correspondingly modified |6| . 

Consider two atoms 1 and 2 at fixed positions (see 
Fig. |l|a), and initially prepared in Stark eigenstates, with 
a dipole moment along z and a given to, as selected by the 
polarization of the laser exciting the Rydberg states from 
the ground state. They interact and evolve according to 
the dipole-dipole potential 
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with r the distance between the atoms. We are interested 
in the limit where the electric field is sufficiently large 
so that the energy splitting between two adjacent Stark 
states is much larger than the dipole-dipole interaction. 
For two atoms in the given initial Stark eigenstate, the 
diagonal terms of Vdip provide an energy shift whereas 
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the non-diagonal terms couple (m, m) ^ (m ± 1, m =F 1) 
adjacent m manifolds with each other. We will assume 
that these transitions are suppressed by an appropriate 
choice of the initial Stark eigenstate Q . As an illustration 
we choose the hydrogen state \r) ~ \n, q ~ n — l,m = 0) 
and find for a fixed distance r — Re^ of the two atoms 
that = (r| (g) (r'|Vdip(^e^)|r) (g) |r) is equal to u{R) = 
—9[n{n — l)]^(ao/i?)'^(e^/87reoao) oc rt^. In alkah atoms 
we have to replace n by the effective quantum number ly 
M. We will use this large energy shift to entangle atoms. 
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FIG. 1. (a) Setup: A constant electric field along the z 
direction is applied to alkali atoms trapped in micro-traps, 
(b) Level scheme: Two ground states \g) and |e) (qubits), 
and laser excitation to the Rydberg state \g} \r). 

The configuration we have in mind is as follows (see 
Fig. 0). We consider two atoms, which for the moment 
are assumed to be at fixed positions Xj with j = 1, 2 la- 
beling the atoms, at a distance R = |xi — X2|. We store 
qubits in two internal atomic ground states (e.g. hyper- 
fine levels) denoted by \g)j = \0)j and \e)j = The 
ground states \g)j are coupled by a laser to a given Stark 
eigenstate \r)j. The internal dynamics is described by a 
model Hamiltonian 



ff*(t,Xi,X2) =w|r}i(r| ® \r)2{r\ 



(2) 



■E 



{S^{t)-t^)\r),{r\ 



^■(i,x,) 



{\9)j{r 



h.c. 



with flj{t,Tij) Rabi frequencies, and Sj(t) detunings of 
the exciting lasers. 7 accounts for loss from the excited 
states \r)j. 

Including the atomic motion, the complete Hamilto- 
nian has the structure 



iJ(i,Xi,X2) = ff^(xi,X2) +i/*(t,Xi,X2) 



(3a) 



the internal Hamiltonian i7'(t,Xi,X2) in Eq. (3b), where 
(the c-numbers) Xj now denote the centers of the initial 
atomic wave functions as determined by the trap, while 
avoiding motional effects arising from H'^{t, Xi, X2). This 
requires that the gate operation time At is short com- 
pared to the typical time of evolution of the external 
degrees of freedom, H'^At <^ 1. Under this condition, 
the initial density operator of the two atoms evolves as 
Pe ® Pi Pe ® p'i during the gate operation. Thus the 
motion described by pe does not become entangled with 
the internal degrees of freedom given by pi. Typically, 
the Hamiltonian will be of the form 
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which is the sum of the kinetic energies of the atoms and 
the trapping potentials for the various internal states. In 
our estimates for the effects of motion we will assume 
that the potentials are harmonic with a frequency lo for 
the ground states, and oj' for the excited state. 

Physically, for the splitting of the Hamiltonian accord- 
ing to Eq. ( |3b| ) to be meaningful we require the initial 
width of the atomic wave function a, as determined by 
the trap, to be much smaller than the mean separation 
between the atoms R. We expand the dipole-dipole in- 
teraction around R, Vdip(f) — u{R) — F{rz — i?) + . . ., 
with F — 3u{R)/R. Here the first term gives the en- 
ergy shift if both atoms are excited to state |r), while 
the second term contributes to H"^ and describes the me- 
chanical force on the atoms due to Vdip. Other contri- 
butions to H'^ arise from the photon kick in the absorp- 
tion \g) — > |r), but these terms can be suppressed in a 
Doppler-free two photon absorption, for example. We 
obtain i/'=(xi,X2) = H^ - i^(f. - R)\r)i{r\ ® |r)2(r|. 

We will now study several models for phase gates 
according to dynamics induced by H^. A schematic 
overview of the internal evolution of the two Rydberg 
atoms is given in Fig. |^ (with shorthand notation \g)i (g 
|e)2 = \ge) etc.). 

Model A: We assume flj ^ u, and in this scheme in- 
dividual addressing of the atoms is not necessary, i.e., 
ill = ^l2 = ^- We set 61 = 62 = 0. We perform the gate 
with three steps: (i) apply a 7r-pulse to the two atoms, 
(ii) wait for a time At = (fi/u, and (iii) apply again a 
TT-pulse to the two atoms. Since the Rabi frequency f2 
is much larger than the interaction energy the first pulse 



i/'^(i, xi, X2) + -H^*(^, xi, X2), (3b) ^ill transfer all the occupation from the states \g)j to 



where describes the motion and trapping of the 
atoms, and Xj are the atomic position operators, and 
define f = Xi — X2 . Our goal is to design a phase gate for 
the internal states with a gate operation time At with 



the states \r)j and the second laser pulse will bring the 
population back to the ground states \g)j. Between the 
two pulses the state \rr) will pick up the extra phase 
(p = uAt. Thus, this scheme realizes a phase gate op- 
erating on the time scale At cx 1/u. We note that the 
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accumulated phase depends on the precise value of u, i.e. 
is sensitive to the atomic distance. The probability of 
loss due to 7 is approximately given by pi = 2lpj/u. Fur- 
thermore, during the gate operation (i.e. when the the 
state \rr) is occupied) there are large mechanical effects 
due to the force F. This motivates the following model. 
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FIG. 2. Schematics of the ideal scheme. The internal state 
\g)j is coupled to the excited state \r)j by the Rabi frequency 
flj{t) with the detuning 5j{t). The state |e)j decouples from 
the evolution of the rest of the system. 

Model B: We assume u 3> ilj. Let us for the moment 
assume that the two atoms can be addressed individu- 
ally, lH], i.e. Vli{t) ^ ^2{t)- We set 5j = and perform 
the gate operation in three steps: (i) We apply a 7r-pulse 
to the first atom, (ii) a 27r-pulse (in terms of the unper- 
turbed states, i.e. it has twice the pulse area of pulse 
applied in (i)) to the second atom, and, finally, (iii) a 
TT-pulse to the first atom. As can be seen from Fig. ||, 
the state |ee) is not affected by the laser pulses. If the 
system is initially in one of the states \ge) or \eg) the 
pulse sequence (i)-(iii) will cause a sign change in the 
wave function. If the system is initially in the state \gg) 
the first pulse will bring the system to the state i\rg), 
the second pulse will be detuned from the state \rr) by 
the interaction strength u, and thus accumulate a small 
phase (f « 7rf22/2u ^ tt. The third pulse returns the sys- 
tem to the state e*^'^^'^-'|.gg), which realizes a phase gate 
with if = TT — (f Pi t: (up to trivial single qubit phases). 
The time needed to perform the gate operation is of the 
order At « 27r/r2i -I- 27r/r22- Loss from the excited states 
\r)j is small provided 'yAt <C 1, i.e. flj ^ 7. If we choose 
u = 1.8GHz, ilj = lOOMHz and 7 = lOOkHz § we find 
a probability of loss from the excited states of pi = 3.4%. 

An adiabatic version of this gate has the advantage 
that individual addressing of the two atoms is not re- 
quired, rii,2(0 = ^{t) and 5i^2{t) = 6{t). In this 
scheme we assume the time variation of the laser pulses 
to be slow on the time scale given by il and 6 (but 
still larger than the trap oscillation frequency), so that 
the system adiabatically follows the dressed states of 
the Hamiltonian W. After adiabatically eliminating 
the state \rr), we find the energy of the dressed level 
adiabatically connected to the initial state \gg) to be 
given by egg{t) = sgn{S) (^S - (5^ + 2rj2)i/2^ /2 with 



6 = S — n"^ /{AS + 2u) the detuning including a Stark 
shift. For the dressed levels connected to \eg) and \ge) we 
have eegit) = sgn{S) {\6\ - {S^ + n^)^/^) /2. The entan- 
glement phase follows as (p{t) = /J dt' {egg{t') — 2€eg{t')). 
For a specific choice of pulse duration and shape Q{t) and 
6{t) we achieve (p = (p{At) — n (see Fig. ||a). In Fig. ||3 
the phases accumulated in the dressed states of \gg) and 
\eg){\ge)), and the resulting entanglement phase are 
shown. To satisfy the adiabatic condition, the gate oper- 
ation time At is approximately one order of magnitude 
longer than in the gate discussed above. 
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FIG. 3. a) The Rabi frequency Q{t) normalized to 
Qo = lOOMHz (solid curve), 5{t) normalized to So = 1.7GHz 
(dashed curve). We chose 7 — lOOkHz, and u = 1.8GHz. b) 
Accumulated phase of the state \gg) (dashed curve) and \ge) 
(dash-dotted curve). The resulting accumulated ip{t) is given 
by the solid curve and the probability of loss from the excited 
state is found to be pi — 1.7%. 

A remarkable feature of model B is that, in the ideal 
limit, the doubly excited state \rr) is never populated. 
Hence, the mechanical effects due to atom-atom interac- 
tion are greatly suppressed. Furthermore, this version of 
the gate is only weakly sensitive to the exact distance 
between the atoms, since the distance-dependent part of 
the entanglement phase (p tt These features al- 

low one to design robust quantum gates with atoms in 
lattices that are not filled regularly. 

We now turn to a discussion of decoherence mecha- 
nisms, which include spontaneous emission, transitions 
induced by black body radiation, ionization of the Ry- 
dberg states due to the trapping or exciting laser fields, 
and motional excitation of the trapped atoms. While 
dipole-dipole interaction increases with i^*, the sponta- 
neous emission and ionization of the Rydberg states by 
optical laser fields decreases proportional to ly'^. For 

< 20 the black body radiation is negligible in compar- 
ison with spontaneous emission, and similar conclusions 
hold for typical ionization rates from the Rydberg states 
for the numbers quoted in the context of Fig. |^. 

We now calculate the motional effects described by i?^ 
on the fidelity of the gate. The dipole-dipole force, given 
by F, causes a momentum kick to both atoms when they 
populate \rr). We assume the atoms to be initially in 
the ground state of the trapping potential. For the adi- 
abatic gate (Model B), we estimate the probability pk 
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of exciting a trap state without changing the internal 
state of the atoms in time dependent perturbation the- 
ory. We find that the perturbative transition probabihty 
is bounded by pk < (3r?flgAt/8u)^ /2 with r/ = a/R < 1. 
For the parameters given in Fig. |[ and = 1/30 we find 
numerically pk ~ 2.5 x 10~^ while the analytic approx- 
imation with a gate operation time At = lOO/rJo yields 
Pk < 2.4 X 10"'^. The probabilities of exciting motional 
quanta, and at the same time changing the internal state 
of the atoms, require the perturbation to induce internal 
transitions with an energy difference of $7 or m in addi- 
tion to causing a motional excitation. Since u ^ ^ w 
whenever there is population in the state \rr) , these prob- 
abilities are much smaller than p^, and are thus negligi- 
ble. For finite temperature of the trapped atoms, we have 
to incoherently sum the probabilities of exciting an atom 
due to the kick for the different trap states. We find 
p^ — {2n -\- l)pk where n is the mean excitation number 
in the trap as determined by the finite temperature T. 

The optical potential seen by the atom in the Rydberg 
state |r) (see Eq. (^) can be different from the trapping 
potential in the ground states \g), |e). This difference 
causes the motional wave function to change its shape 
while the atom is in |r), and thus excites the atomic 
motion. We treat the deviation of the potential in the 
excited state from the ground state trapping potential as 
a perturbation and estimate the probability pt of excit- 
ing an atom from the vibrational ground state during the 
gate operation. We find that the perturbation theory ex- 
pression is bounded by pt < — ci;'^|At^/128. For the 
parameters given in Fig. ^, and trap frequencies of w = 
IMHz and oj' = 500kHz we find numerically pt ~ 10~^. 
The analytical approximation (with At = lOO/rio) gives 
Pt < 3.9 X 10~^. At finite temperature T we find by 
incoherent summation = {2n? + 2n -\- l)pt- 

Alternatively, the optical trapping potential can be 
turned off for the short time of the gate operation. The 
shape of the wave function of the atoms evolves then in- 
dependently of their internal states |e}, \g) or \r). There- 
fore, no entanglement between external and internal de- 
grees of freedom is created during the gate operation. 
However, the releasing and retrapping will cause heating 
of the atoms. We estimate this effect at finite tempera- 
ture T and find an increase in the mean excitation num- 
ber of n ^ h + An — n + (wAt)^ (2n-|- 1)/4 during a gate 
operation. 

The influence of imperfections on the other schemes 
discussed in this Letter can be estimated in the same way 
as for the adiabatic gate. However, in model A the mo- 
mentum kick will be stronger than for the other schemes 
because all of the population is transferred to \rr) for a 
time if/u. The force experienced by the atoms during 
this time is approximately given by F, yielding a total 
momentum transfer of 3ip/R to the atoms, which for a 
separation i? ^ A is of the order of several recoils hk with 



k = 2tt/X. Therefore, only if the atoms are confined deep 
in the Lamb-Dicke regime 27ra/A <C 1 with A the opti- 
cal wave length, will this momentum transfer not cause 
significant excitation. 
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[10] In fact, this dependence can be entirely eliminated in the 
approach where splitting rather than shift of the state 
\rr) is involved. This regime can be achieved by choosing 
DC electric fields so that the energy splitting of a pair of 
transitions np —> ns and np {n + l)s becomes degen- 
erate (see e.g. I. Mourachko, et al., Phys. Rev. Lett. 80 
253 (1998); W.R. Anderson et al., ibid., 249), which re- 
sults in rapid hopping of excitation between atoms and 
associated energy splitting. See also M.D. Lukin and 
P.R. Hemme, Phys. Rev. Lett. 84, 2812 (2000). 
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